Phase estimation as a quantum nondemolition measurement 
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The phase estimation algorithm, which is at the heart of a variety of quantum algorithms, includ- 
ing Shor's factoring algorithm, allows a quantum computer to accurately determine an eigenvalue of 
an unitary operator. Quantum nondemolition measurements are a quantum mechanical procedure, 
used to overcome the standard quantum limit when measuring an observable. We show that the 
phase estimation algorithm, in both the discrete and continuous variable setting, can be viewed as 
a quantum nondemolition measurement. 
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There is currently a great deal of active research aimed 
at developing implementation schemes for quantum com- 
putation. For an overview of this field see Nielsen and 
Chuang [[tj or Preskill notes and references therein. 
The goal of quantum computation is to build a com- 
puter, which, using the laws of quantum mechanics, can 
out perform any classical computer. Using the laws of 
quantum mechanics to enhance a device's performance is 
not new. In quantum optics, the procedure of quantum 
nondemolition (QND) measurement || ||, ||, ||[ [?J aims to 
accurately measure an observable without perturbing it. 
This procedure is consistent with the uncertainty princi- 
ple, as, ideally, noise is induced only in the conjugate of 
the observable being measured. 

Much of the impetus for research in quantum com- 
putation stems from Shor's algorithm which allows a 
quantum computer to factor integers exponentially faster 
than can currently be done on a classical computer. The 
factoring of large integers allows the decryption of many 
public key encryption systems. It has been shown that 
Shor's, and many other quantum algorithms, can be de- 
scribed in terms of the phase estimation algorithm ||, [l0| . 

It is evident that both phase estimation and QND mea- 
surement have the same goal: they both aim to accu- 
rately measure the eigenvalue of an operator (observ- 
able). In this paper we show that they are essentially 
the same procedure. That is, every instance of the phase 
estimation algorithm can be considered a QND measure- 
ment. 

The original motivation for the study of QND measure- 
ments was the need to accurately measure weak classical 
forces ||. Such measurements will be necessary, for ex- 
ample, in the detection of gravitational radiation. Clas- 
sically, it is possible to measure both the position and 
momentum of a particle with complete precision. How- 
ever, quantum mechanics tells us that the product of the 
uncertainty in the position and momentum of a particle 
must be greater than Planck's constant divided by 4-7T, 



cillator is the fact that any quantum mechanical state is 
characterized by an "error box" in the phase plane with 
area at least 



AxAp/muj > h/2mtu, 



(2) 



where m and to are the mass and angular frequency re- 
spectively [||. We can re- write Eq. (g) in terms of the 
quadrature variables, X and Y, 



where 



AX AY > h/2r 



+ ip/mtu = (X + iY)e 



(3) 



(4) 



If we naively perform an amplitude- and-phase measure- 
ment of the quadrature variables, we obtain a round error 
box, with 



AX = AY = 



2mu) 



(5) 



However, by performing a back-action-evading QND 
measurement of, say, X, the error box becomes a long 
thin ellipse, with 



AX < AF 



(6) 



and AX can be made as small as one wishes, in principle. 

Although a precise definition for a QND measurement 
is hard to pin down from the literature, there are three 
basic criteria which a QND measurement must satisfy. 
Firstly, it is required that the QND variable A, which we 
wish to measure, commutes with the free Hamiltonian 
evolution of the system, 



[Hp, A] = 0, 



(QND 1) 



AxAp > -. 



(1) 



so that if the system is measured in an cigenstate of A, it 
remains in this cigenstate for all subsequent times, (per- 
haps only in the interaction picture). The most impor- 
tant criterion which must be satisfied, is that the Hamil- 
tonian, Hsm, coupling the meter variable B to the sys- 
tem must commute with the QND variable A, 



The consequence of the uncertainty principle for an os- 



[Hsm,A]=0. 



(QND 2) 
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This is known as the back action evasion criterion. Fi- 
nally, it is also necessary that the meter variable B, does 
not commute with the interaction Hamiltonian, 



[Hsm,B]^Q. 



(QND 3) 



If this requirement is not fullfillcd. then measuring the 
meter variable B will yield no information about the 
QND variable A. 

To illustrate how a QND measurement procedure 
works, consider the example of an ideal quadrature phase 
measurement of a single mode field 11 . Suppose we 



have two degenerate modes a and b, of an electromagnetic 
field. We wish to measure X a using the meter variable 
Xj). Now, if we assume that the amplitude quadratures 
of each mode can be coupled according to the interaction 
Hamiltonian, 



H. 



SM 



hxX a Y b . 



(7) 



where x is the coupling strength, then clearly (QND 2) 
and (QND 3) arc satisfied. (QND 1) is also satisfied by 
choosing a rotating frame of reference. In order to draw 
analogies with the phase estimation algorithm, we will 
divide the QND measurement into three steps. The first 
step involves initializing mode b into the zero eigenstate 
of the meter variable, \0)b, 



X b \0) t 



0. 



(8) 



Now, for simplicity, we will assume that mode a is already 
in an eigenstate \x) a of the QND variable, 

X a \x) a = x\x) a . (9) 

Thus, the initial state of the quantum system is 

\x) a \0) b . (10) 

Re- writing the meter system in terms of the phase vari- 
able, Yb, gives 



dy \x)a\y)b- 



(ii) 



Evolving according to the interaction Hamiltonian 
(Eq. (^)) results in the state 

e ixX u Y b fdy\x) a \y)b = f dy e?™\x) a \y) b . (12) 



Finally re- writing mode b in terms of the amplitude vari- 
able Xb, we have 



\x) a \xx)b- 



(13) 



Therefore, by performing a quadrature measurement of 
the amplitude of mode b, we obtain information about 
the state of mode a. 



In the above description we have used eigenstates of the 
quadrature variables, however, like position and momen- 
tum eigenstates, these are physically unrealisable. How- 
ever, they can be approximated to arbitrary degree by 
using highly squeezed states. 

The goal of the phase estimation algorithm is to obtain 
an eigenvalue of a unitary operator U . Although hybrid 
systems have been discussed in the literature [jl2[ [l3| , 
generally the phase estimation algorithm is described in 
terms of (discrete) qubit registers |l(| . Discrete classi- 
cal computation, is accomplished through the manipula- 
tion of binary digits (bits) using logic gates such as AND, 
OR, NOT and FANOUT, whereas discrete quantum com- 
putation uses arrays of two level quantum systems (qubit 
registers) which are manipulated via one and two qubit 
unitary operators. In analogy with the classical bit, it is 
conventional to use the symbols |0) and |1) to denote the 
two levels of the quantum system. Also in analogy with 
classical computation, we use the compact integer repre- 
sentation to represent the basis states of a qubit register. 
For example, 



|l)s®|l)a®|0)i®|l)o = |13). 



(14) 



For an m qubit register, we also introduce the Hermitian 
operator, 



m— 1 

A = ^2 I ® ■ ■ ■ ® I ®2 i ~ 1 (I - Z)®I < 

< 7 tn term 



■ 1(15) 



where / is the identity operator, and Z is the pauli op- 
erator, with matrix representation, 



1 

-1 



(16) 



The computational basis states are eigenstates of the op- 
erator A with an eigenvalue corresponding to their integer 
representation, 

A\x) = x\x) Vi6 [0,l,...,2 m - 1]. (17) 

Another basis which is useful in understanding the phase 
estimation algorithm is the Fourier transformed (FT) ba- 
sis. This is the basis which arises by performing the quan- 
tum Fourier tranform fl4jl , Q, on the computational basis, 

\x) = Q\x) Vi e [0,1,..., 2 m - 1]. (18) 

The Hermitian operator associated with this basis is 



as 



QAQ~ 



T\x) = x\x) Vie [0,1, 



(19) 



(20) 



Suppose that we have a register of qubits, the target 
register, which is in an eigenstate \4>)t, of some unitary 



3 



operator U . We wish to determine the eigenvalue corre- 
sponding to \4>)t- As U is unitary, we have 



- e ir t>\ 



(21) 



where <p is some number between and 2-7T. The phase es- 
timation algorithm utilizes an auxiliary register of qubits, 
which we shall call the index register. Fig. |l| depicts the 
quantum circuit diagram for the discrete phase estima- 
tion algorithm. The algorithm begins by placing the in- 



" \o>, -m- 



\o> -m — ^ 




FIG. 1: Quantum circuit diagram representing the discrete 
phase estimation algorithm. 



dex register in the zero of the FT basis, 
|*o) = |O)/|0)t, 



(22) 



Assuming that each qubit in the index register begins the 
computation in the state |0), this step is accomplished 
by performing a Hadamard gate, H 7 on each qubit. The 
next step is to perform a series of controlled U 2 gates, 
which couple the target and index registers. For each 
k e [0,1,..., vn — 1] , the gate U 2 is applied condition- 
ally upon the state of the A:-th qubit in the index register. 
The cummulativc effect of applying these gates is the op- 
eration, U A . The state of the system after the application 
of U k is 



^ A |*o> 
1> A |O),|0) T 
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(23) 



3=0 



where the last line makes use of Eq. (^l|). The final step 
in the algorithm involves measuring the index register in 
the FT basis, or equivalently, performing the operation 
and measuring in the computational basis. This 
will, with high probability G3|, result in an m-bit ap- 
proximation of </>. 

Having re-cast the phase estimation algorithm in terms 
of the FT basis, it is now not difficult to see that the phase 



estimation algorithm is effectively a QND measurement. 
We are trying to find an eigenvalue of the unitary oper- 
ator, U, which can of course be written as 



U 



(24) 



for some Hermitian operator, Hjj. This operator is the 
QND variable that we are trying to measure, T is the 
meter variable, and HjjA is the interaction Hamiltonian, 
as 



AHrrA 



(25) 



In the quantum circuit model, it is assumed that the 
qubits undergo no free evolution, so (QND 1) is auto- 
matically satisfied. It is clear from Eq. (||^) that the 
back action evasion criterion is satisfied (QND 2), and 
it is not difficult to see that the interaction hamiltonian 
does not commute with the meter variable (QND 3) as 



[A,f] 0. 



(26) 



To illustrate how the phase estimation algorithm acts as 
a QND measurement, we describe the simplest possible 
algorithm. Suppose we wish to determine the eigenvalue 
of the pauli X operator, or NOT gate, 



X = 



1 

1 



(27) 



The eigenstates of this operator (ignoring a normalization 
factor) are |0)±|1). The quantum circuit diagram for this 
simple phase estimation algorithm is depicted in Fig. ||. 
In this example the QND variable is 



I l0> " 



a \o> ± \i> 
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\o> ± li> 



FIG. 2: A simple example of the phase estimation algorithm 



Hu — —{It — X T ) 1 
the meter variable is 

f = QAQ' 1 = 1 



1 -1 
-1 1 



(28) 



(29) 



as the single qubit quantum Fourier transform is the 
self inverse Hadamard operator, Hi, and the interaction 
Hamiltonian is 



H V A = -(I T - X T ){h - Zi). 



(30) 
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Obviously, all the QND criteria are satisfied, and the 
measurement proceeds by placing the meter system into 
the zero eigenstate of the meter variable, |0) + |1), which 
is accomplished by the first Hadamard gate, evolving ac- 
cording to the interaction Hamiltonian, which is accom- 
plished by the controlled NOT gate, and finally measur- 
ing the meter variable, by performing another Hadamard 
and then measuring in the computational basis. 

We have just shown that the discrete phase estima- 
tion algorithm is equivalent to a QND measurement. 
Recently, a number of papers have been written ex- 
ploring quantum computation using continuous variables 
111, 111, ~ 



18 1 . We now show that continuous vari- 



able phase estimation is equivalent to continuous variable 
QND measurement. 

As with all quantum algorithms, the phase estimation 
algorithm can be partitioned into three stages; the initial- 
ization stage, the entangling stage and the measurement 
stage. The initialization stage involves placing the in- 
dex system into the zero eigenstate of the FT basis. The 
entangling stage involves the application of the operator 



iHuA 



and the measurement stage involves measuring the index 
system in the conjugate of the computational basis. If 
we replace the discrete qubit registers with infinite level 
(continuous) quantum variables, and carry out each of 
these steps, we obtain something very similar to the ideal 
quadrature QND measurement discussed earlier. Sup- 
pose that the continuous variable operator, acting on the 
target system, whose eigenvalue we wish to determine 
is e lAT , and we choose the position of the continuous 
index system as the computational basis. Then follow- 
ing the steps of the phase estimation algorithm, we pre- 
pare our index system in the zero eigenstate of the conju- 
gate basis, that is, the zero momentum eigenstate. This 
step corresponds to initialising our meter variable. Then 
we perform the U operator, conditioned upon the com- 
putational basis, e lXlAT , which corresponds to applying 
an interaction Hamiltonian which satisfies (QND 1) and 
(QND 2). Before finally performing a Fourier transform 
and measuring the position variable, which is equivalent 
to measuring the momentum of the index system (i.e. the 
meter variable). This is clearly a QND measurement of 
A. 

Of course, position eigenstates are physically unrealis- 
able, and any continuous variable implementation of the 
phase estimation algorithm would need to be done us- 
ing approximate position eigenstates, such as squeezed 
states. Just as ideal QND measurements can only be 
approximated experimentally, ideal continuous variable 
phase estimation can also only be approximated. 



The study of quantum computation brings together a 
number of disciplines including computer science, mathe- 
matics and quantum mechanics. We have taken the phase 
estimation algorithm, which has its roots in computer 
science, and shown that it is equivalent to the quan- 
tum mechanical procedure of QND measurement. We 
have also described a continuous variable version of the 
phase estimation algorithm. It is interesting to speculate 
on whether other algorithms such as Grover's algorithm 
might have similarities with well studied quantum phe- 
nomenon. 

BCT and GJM thank S. Lloyd for helpful discussions. 
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